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ABSTRACT 



Herschel's PACS instrument observed the environment of the binary system Mira Ceti in the 70 and 160 jim bands. These images 
reveal bright structures shaped as five broken arcs and fainter filaments in the ejected material of Mira's primary star, the famous 
AGB star o Ceti. The overall shape of the IR emission around Mira deviates significantly from the expected alignment with Mira's 
exceptionally high space velocity. The observed broken arcs are neither connected to each other nor are they of a circular shape; they 
stretch over angular ranges of 80 to 100 degrees. By comparing Herschel and GALEX data, we found evidence for the disruption of 
the IR arcs by the fast outflow visible in both Ha and the far UV. Radial intensity profiles are derived, which place the arcs at distances 
of 6-85" (550 - 8000 AU) from the binary. Mira's IR environment appears to be shaped by the complex interaction of Mira's wind 
with its companion, the bipolar jet, and the ISM. 
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1. Introduction 

Centuries ago o Ceti confused the astronomers all over the world 
with its appearance and disappearance o n the nigh t sky a nd was 
therefore named Mira, the wonderful, by lHeveliusl (Il662l) . 

Today we know that Mira is a wind-accreting binary star sys- 
tem (iReimers & Cassatellalll985l) at a distan ce of 91.7 pc (cor- 
responding to a parallax of 10.9 ±1.2 mas: Ivan Leeuwenll2007l) 
consisting of the primary star Mira A (o Ceti) and its compan- 
ion Mira B (VZ Ceti). Mira A is an M-type oxygen-rich star 
with a pulsation period of 331 d and a brightness amplitude in 
the visible of 8 mag. Moreover, Mira A is the prototype of stars 
in their late stage of evolution residing on the asymptotic gi- 
ant branch (AGB). The nature of Mira B has been unknown for 
long; with an effective temperature of about 100 00 K, it might 
either be a main- sequence star or a white dwarf (iKarovska et al.l 
11997b . but recent observations show strong evidence in favour 
of a white dwarf (iSokolos ki & Bildsten 201 (J) .From a Hubble 
Space Telescope observation, IKarovska et al.l (1 19971) found the 
separation of Mira AB to be ^0.6", which corresponds to 55 AU. 
A very uncertain orbit, based on astrometric measurements since 
1923, which has a period of ^500 yr, has been published by 
iPrieur etal.1 (120021) . 

Stars on the AGB loose a significant fraction of their mass in the 
form of a stellar wind. For Mira A, a moderate mass loss rate of 
10" 7 M© yr" 1 was determined by Mauron &Cauxl (119921) . which 
corresponds to a dust mass-loss rate of about 5.5 x 10" 10 M© yr 



(Meixner, 201 1, priv. comm.). From mm-CO lines Jjosselin et al.l 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



(2000) found a wind expansion velocity of 4 km s 1 and a 
molecular envelope around the binary system with a size of 
about 20". This envelo pe is disrupted by a slow bipolar out- 
flow of about 8 km s" 1 dJosselin et al.ll2000l) . but a high- velocity 
com ponent at 160 km s" 1 has been observed by Mea burn et al.l 
( 2009). In the outer environment of Mira, a spectacular obser- 
vation with the GALEX ultraviolet satellite dMartin et al.ll2007l) 
revealed a cometary-like structure wit h an extension of m ore 
than 2° (« 4 pc) from head to tail. iMartin et all (120071) at- 
tribute this tail to Mira's exceptionally high space velocity of 
^ 1 10 km s" 1 and the resulting strong interaction of the stellar 
wind with the interstellar medi um (ISM). This structu re was also 
observed in the 21 cm H I line (Matt hews et al . 2008) a nd is well 
reproduced by the hydrodyna mic simulations of IWareing et al.l 
(120071) . iRaga & Cantdl (120081) and lEsquivel et al] (I201QI) . In the 
vicinity of the binary system, knot-shaped enhancements of the 
UV and Ha emission were found at a distance of V to 2' from the 
binary dMartin et al.ll2007l; iMeaburn et al.ll2009l) . corresponding 
to 5500 to 11 000 AU, and in a direction cons i stent w ith that of 
the bipolar outfl ow observed by lJosselin et al.l (l2000h . In the far 
IR lUetal (120081) resolved with the Spitzer Space Telescope the 
contours of Mira's astropause including the circumbinary enve- 
lope, termination shock, and astrosheath. 

The present study offers a detailed view of Mira's circumstel- 
lar shell(s) thanks to the high spatial resolution of the PACS in- 
strument (which has a point- spread function of 576 full- width 
at half-maximum at 70 fim) onboard the Herschel satellite. The 
images presented here reveal clumpy broken arcs (see Sect. 13.11) 
that result from the complex interactions of the wind with the 
ISM, the companion and the fast outflow. 
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Fig. 1. Panel (a): deconvolved PACS image at 70 jim. The ar- 
row indicates the space motion and the position in 500 yrs; (b) 
is the same for the 160 jim band and (c) is the 70 yum decon- 
volved PACS image with contours and arcs labelled as A, B, 
C, D and E. The lowest contour is the 3<x threshold (and is de- 
picted in white), the arrow labelled with V space shows the di- 
rection of the space motion, the bar labelled PA orbit gives the 
orientation of the major axis of the apparent orbit, orthogonal to 
the node line. Panel (d) results from the 'toy model' described 
in Sect. 13.31 that is based on the hvdrody namical simulations of 
iMohamed & Podsiadlowskil d20Q7ll20Tlb , 

Table 1. Kinematic data for Mira. Negative inclinations are 
towards the observer. The various vectors are displayed in 
a 3D view in Fig. [4] The adopted solar motion to con- 
vert heliocentric velocities in LSR velocities is (U, V, W) Q = 
(11.10, 12.24,7.25) km s" 1 (iSchonrich et alJl2Q10b , The matrix 
to co nvert equatorial velocities int o galactic velocities is taken 
from lJohnson & Soderbloml (119871) . 



b(°) 


-58 


tjj (mas) 


10.9 + 1.2° 


Va,i£R (mas yr ') 


-20.93° 
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-222.39° 
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a: vanLeeuwen (2007); b: Mauron & Caux (1992); c: Josselin et al. 
(20001): d : iMartin et all 00071) : iMeaburn et all d2009): e: iPrieur et al.l 
(2002); f: lFong et all 1*2006) 



2. Observations and data reduction 

The observations of Mira took place on February 9, 2010, 
using the PhotoArr ay Camera and Spectrometer (PACS; 
iPoglitsch et~ai1 l2010h onboard Herschel as part of the MESS 
(Mass-loss of Evolved S tarS) Guaranteed Time Key Program 
dGroenewegen et al]|201 ll) . 

The PACS camera observed at 70 yum and 160 yum 
simultaneously, covering a field of 30' x 30' on the 
sky with a scan-map mode. For the basic data pro- 



cessing to produce the images displayed in Fig. [T] the 
Herschel Interactive Processing Environment (HIPE) was used 
(lArdila & Science Ground Segment Consortium! 12010) together 
with Scanamorphos, an IDL software for building sc an maps 
from PACS (and SPIRE) observations (lRousselll20lTh. Finally 
we d econvolved the images as described by lOttensamer et al.l 
(1201 ll) . Although the PACS instrument offers a resolution of 
3.2" per pixel in the 70 jim (blue) band and 6.4" per pixel at 
160 yum (red band), the images shown in Fig. [T] are oversampled 
by a factor 3.2, resulting in a sampling of 1" and 2" per pixel (at 
70 and 160 yum, respectively), to be compared with the FWHM 
of 576 and 1 174 at those wavelengths. We estimate the average 
background emission to be 1.5 x 10" 4 Jy pixel" 1 in the blue and 
3 x 10" 4 Jy pixel" 1 in the red band, yielding the l<x sensitivity. 

3. Results and discussion 

3.1. I R images 

Mira's environment is remarkably structured and exhibits arc- 
like structures in the surrounding of the binary system, labelled 
A to E in Fig.[T]:. The best-defined and brightest arc, C, stretches 
over position angles (PAs) ~ 190° -270° (south-west to west), is 
well detached from the star and covers a range of angular separa- 
tion from 67" (6100 AU) at the closest point to 85" (7800 AU) 
at the most distant point from the star. The intensity was inte- 
grated over this PA range and the resulting radial intensity pro- 
file (Fig. [2]) clearly shows the arc as a plateau in the diagram. The 
extension of this plateau confirms that the arc is non-circular or, 
if circular, is not centred on the position of the star. We found an- 
other arc, D, in the same direction as arc C but much closer to the 
star and stretching over much smaller PAs. In the radial intensity 
profile, arc D is visible as a bump located at a distance of 21- 
28" (1900-2600 AU). East of Mira more arcs are visible. In the 
south-east direction, we see arc B close to the star. It stretches 
over PAs 75° to 180° and is more circular than arc C and there- 
fore appears as a narrow peak in the radial intensity profile at a 
distance of 41" (3800 AU) from the binary (Fig.0. 
Farther out, directly to the east, we find another structure, arc A, 
stretching towards north-east. In the integrated intensity profile, 
it can be recognized as a small plateau representing its fairly faint 
nature at a distance of 68-89" (6200-8200 AU). This matches 
the distance of arc C, which leads to the suggestion that both arcs 
are somehow related to each other. 

It is also possible that another arc, E, is present in the south-east 
direction, very close to the binary at a distance of less than 15" 
(1400 AU), of which we only see the edge at a PA of about 65°. It 
is noteworthy that the arc sequence E-D-B-C-A is a sequence of 
both increasing distance from the star and increasing PA range, 
arc E being the shortest and closest to the star, and A the longest 
and farthest. 

The overall 3cr emission above the 70 yum IR background of 
Mira, visible in FigQ]: as the outermost contour, spans as far 
as 157" towards the south-east. It is remarkable that the exten- 
sion of the emission is towards a position angle of « 160° be- 
cause one would expect the emission to be aligned with Mira's 
space motion, but that motion is in the direction of PA = 185°. 
In Fig. [T]: we also indicated the orientation of the major axis 
of the apparent or bit, which is oriented along a PA of ^ 1 40° 
(iPrieur et al.l 120021) . As discussed in more detail in Sect. 13.31 
the orbital motion and the mass transfer from the binary can 
also shape its environment. Nevertheless, the extension of the 
IR emission lies between both values and makes a distinct attri- 
bution difficult. 
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Fig. 2. Radial intensity profiles of Mira in the 70 yum band as a 
function of the distance to the star. Continuous (blue) line: in- 
tegrated intensity across the arcs B and A (75° < PA < 180°). 
Dashed (red) line: across the arcs D and C (190° < PA < 270°). 
The top scale gives the distance from the star computed from 
Mira's distance of 91.7 pc. 

3.2. Comparison with GALEX UV observations 

Before drawing conclusions from the structures observed in 
the IR, it is interesting to confront these IR observations with 
GALEX ultraviolet observations of Mira in the far and near UV 
(152.8 nm and 2 27.1 nm, respectivel y; iMartin et all l2007h . In 
the far UV image, Marti netaD d2007) found emission knots ar- 
ranged north-west and south (Fig. 0). The axis of the streams 
is consistent with the direction of the bipolar outflow detected 
by iJosselin et al.1 (120001) . and moreover, the U V knots are ex- 
actly coincident with similar knots seen in H a (Mea burn et al.l 
12009b . In contrast, Fig. [3] reveals that in the PACS 70 jim im- 
age the southern stream exactly corresponds to a region with an 
IR flux below the 3<x thresholcQ. It thus looks as if the south- 
ern stream is penetrating and disrupti ng the material resp onsi- 
ble for the IR emission. According to Mea burn et al.l (120091) . the 
southern stream is indeed inclined at an angle of 69° with re- 
spect to the plane of sky towards us. Put together, all these ar- 
guments clearly indicate that the IR emission must come from a 
3D structure and is not confined to the orbital plane; otherwise, 
there could be no explanation for the anti-coincidence between 
the GALEX far UV/Ha emission and the IR emission. 
Another interesting clue is that the northern lobe of the bipo- 
lar flow remains visible even though it is located behind the IR 
shell. The dust present in the shell will efficiently absorb the far 
UV radiation, so that the intensity of the northern lobe measured 
by GALEX should be smaller than that of the southern lobe, 
which has no intervening dust. In the far UV, the peak intensity 
of the northern lobe is indeed about 1.5 times fainter than the 
peak intensity of the southern lobe. Assuming that the northern 
and southern jets carry the same energy, this ratio may then be 
seen as resulting from dust extinction: 

1/1.5 = exp(-W£ext,uv), (1) 

where N is the dust column density, a is the grain radius, 
and <2ext,uv is the ratio between the extinction cross section 

1 The 'hole' in the IR emission located about V to the East is less 
significant than the one considered, with fluxes about twice smaller for 
the latter. 




Fig. 3. GALEX far UV image of Mira obtained bv lMartin et al.l 
(120071) overplotted with contours of the Herschel 70 yum image, 
the thick blue contour corresponding to an IR flux just 3cr above 
background. The bright white regions correspond to strong UV 
emission at both ends of the fast jet. Note the absence of IR 
emission at the end of the southern jet. The red arrow indicates 
the position of the star in 500 yrs. The bow shock is the faint 
structure in the background in the south-west. 

and the geometrical section. This constraint can be inserted in 
Kirchhoff 's law for the IR emission j IR , leading to 

JlR = #IR(7dust) ln(1.5) Qabs,IR/Gext,UV, (2) 

where #ir is the Planck intensity function at the dust tempera- 
ture and <2abs,iR the ratio o f the absorptio n cross section and the 
geometrical section. From iDrain el (1 1985b . g a bs,7(W2ext,o.i5/;m ~ 
10" 3 (almost independent of the grain size), which leads to 

j IR = 4.2xlO- 4 5 IR (r dust ). (3) 

From the measured intensity ratio at 70 and 160 jum for the 
northern lobe, we find a dust temperature of 36 K; combining 
^70//m(36 K) with the measured jVo^m = 4.42 x 10 7 Jy/sr, we ob- 
tain j70//m/#70//m(36 K) = 1.2 x 10" 4 , iwhich agrees well with 
the value predicted by Eq. [3] 

3.3. Roie of the binary system 

The arcs observed aro und Mira (Sect. 13.11) are (a lmost) unique 
in the MESS sampl<£| dGroenewegen et aljfcoi lL Cox et al., in 
preparation), which makes it tempting to relate them to the bi- 
nary nature of Mira. 

In binary systems involving a mass-losing component, 
nested half- shell s are seen when the orbit is orientated close 
to edge-on dHel 120071) . and they represent the 3D structure 
extending above the Archimedian spiral present in the plane 
of a circular orbit (for eccentric orbits, the spiral becomes 
broken). T he latter was predic t ed by t he hydrodynamics sim- 
ulations of [j heuns & Jorissen (1993), Mast rodemos & Morris 
(1999), Mohamed & Podsiadlowski ( 2011), and detected around 

2 The prototypical carbon star IRC +10216 is surrounded by a large 
num ber of thin and faint arcs, known for quite some time in the op- 
tical (Mauron & Huggins 1999, 2000), and now identified in Herschel 
images as well (Decin et al., in preparation). Their appearance is, how- 
ever, quite different from the few thick arcs seen around Mira, in line 
with the fact that t here is so far no ev idence that IRC +10216 belongs 
to a binary system (Huggins et al. 20091). 
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Fig. 4. Various vectors entering the discussion in a 3D view, with 
the red and green planes corresponding to the two possible or- 
bital orientations. Because bipolar flows are generally perpen- 
dicular to the orbital plane, the green orientation seems the most 
likely. The sky plane is depicted in blue. The thick arrows are 
projections onto the sky plane. The observer is located on the 
bottom plane and looking upwards. 

AFGL 3068 by iMauron & Hugginsl (120061) and iMorris et al.l 
(2006). In Mira, the observed arc structure resembles neither 
an Archimedian spiral, nor nested half-shells; its interpretation 
must therefore be more complex, involving two supplementary 
physical processes at play, namely the presence of a bipolar jet, 
and the interaction with the ISM (to fix the ideas, Fig.|4]presents 
the respective orientations of all relevant vectors involved in the 
Mira system). We already arg ued in Sect. 13. 21 that the fast bipo- 
lar outflow (160 km s" 1 : lMartin et al .120071 iMeaburn et al.l2009l) 
seems to cut a hole in the IR shell (see Fig. [3]). Similarly, it could 
easily disrupt the arcs with its high velocity, and may actually be 
responsible for cutting the edges of arcs B, C, and D, because all 
these edges are located along the southern jet, at a position angle 
close to 180°. 

On the other hand, it is likely that the interaction with 
the ISM, more precisely, the termination shock, plays some 
role in shaping (at least) arc C, as already suggested by Ueta 
(2008) from lower-resolution Spitzer images, and confirmed in 
the Herschel images by the sharp ridge seen on the contour lev- 
els in Figs. [J: and [3] Arc C is indeed located almost upstream 
along Mira's space motion. On the 160 jim image (Fig. [TJ)), 
arcs A and C extend towards the north and join with each 
other, exactly opposite to the space motion of Mira, in a tail- 
like or bull et-like structure, as seen in the hydrodynamic simu- 
lations bv lWareing et al.l (12007b and lEsquivel et al I (l2Q10h . This 
is clear evidence that the outer arcs feel the space motion and 
interact with the ISM. Another indication thereof is provided 
by the curva ture of the bipolar-outflow southern stream, which 
Mar tin et all (120071) attribute to its deceleration by the post- 
shock flow. Therefore, the arcs seen around Mira likely result 
from a combination of the projected 3D structures resulting from 
the interaction of Mira's wind with its companion on one hand, 
and with the ISM on the other hand. 

As a fi rst exploration of this combination , we us ed the sim- 
ulations of lMohamed & Podsiadlowskil (I2007LI201 lb of an AGB 
star nearly filling its Roche lobe in a system 10 AU-wid^E and 
losing mass through a slow wind. The 3D structure resulting 
from the expanding wind disturbed by the companion's gravi- 

3 Although this separation is much smaller than the actual orbital 
separation of Mira system, this difference should not impact our con- 
clusions, because the spiral properties are largely scale-invariant. 



tational pull was then intersected with a paraboloid represent- 
ing the bow shock. All particles within a thin shell along the 
paraboloid surface were kept and the resulting particle density 
projected onto the sky plane (the respective orientations of the 
paraboloid, the orbital plane, the sky plane, and the compan- 
ion were taken into account, as represented in Fig. |4]). If the IR 
emission observed by Herschel around Mira results from the re- 
heating of the wind by the shock, this 'toy model', presented in 
Fig. [T]l, should then give a fair account of the situation. Fig. QJi 
shows that there are indeed promising resemblances between the 
Herschel image and the one from the 'toy model' , apart from the 
absence of brightening of the south-west arc C (where the com- 
pression from the ISM - not included in our model - is maxi- 
mum), and apart from the ISM- swept appearance of the outer- 
most northern arc. Hydrodynamical simulations of mass transfer 
in the Mira system, including the ISM interaction, are thus badly 
needed. 
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